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Fibrillar aggregates of misfolded amy- 
loid proteins are involved in a variety 
of diseases such as Alzheimer disease 
(AD), type 2 diabetes, Parkinson, 
Huntington and prion-related diseases. 
In the case of AD amyloid (3 (A(3) pep- 
tides, the toxicity of amyloid oligo- 
mers and larger fibrillar aggregates is 
related to perturbing the biological 
function of the adjacent cellular mem- 
brane. We used atomistic molecular 
dynamics (MD) simulations of A|3 , g 
fibrillar oligomers modeled as protofila- 
ment segments, including lipid bilayers 
and explicit water molecules, to probe 
the first steps in the mechanism of 
A|3-membrane interactions. Our study 
identifies the electrostatic interaction 
between charged peptide residues and 
the lipid headgroups as the principal 
driving force that can modulate the 
further penetration of the C-termini 
of amyloid fibrils or fibrillar oligomers 
into the hydrophobic region of lipid 
membranes. These findings advance our 
understanding of the detailed molecular 
mechanisms and the effects related to 
A|3-membrane interactions, and suggest 
a polymorphic structural character of 
amyloid ion channels embedded in lipid 
bilayers. While inter-peptide hydro- 
gen bonds leading to the formation of 
|3-strands may still play a stabilizing 
role in amyloid channel structures, these 
may also present a significant helical 
content in peptide regions (e.g., termini) 
that are subject to direct interactions 
with lipids rather than with neighbor- 
ing A|3 peptides. 



Introduction 

Alzheimer disease (AD) is characterized by 
a strong loss of cognitive functions that is 
pathologically correlated with the appear- 
ance in neural tissue of fibrillar deposits 
in extracellular plaques and intracellular 
tangles containing amyloid P (Ap) pep- 
tides and tau proteins, respectively. Two 
mutually non-exclusive concepts emerged 
regarding Ap neurotoxicity, one is that the 
peptides may be cytotoxic by themselves 1,2 
and that their aggregation is required for 
toxicity, 3,4 and the other is that the amy- 
loid peptides may interfere synergistically 
with the normal biochemical processes 
involving other molecules. 5,6 

It has been demonstrated that amyloid 
plaques can be induced in mice via intra- 
cerebral 7 or intra-peritoneal 8 injection of 
misfolded AP-containing brain extract, 
this behavior of the peptide being similar 
to that of prions. 9 However, there is no 
direct evidence that injected Ap is both 
necessary and sufficient to trigger cere- 
bral amyloid deposition without cofactors 
that facilitate its aggregation. Induction of 
amyloid formation in vivo with synthetic 
Ap has been unsuccessful so far, leading 
to the hypothesis that an amyloid-enhanc- 
ing factor or a particular peptide confor- 
mation is required to trigger its deposition 
in vivo. 7 

Several studies have supported the 
idea that the interactions between Ap 
peptides and cellular membranes lead to 
the formation of fibrillar Ap aggregates 
and to their subsequent cytotoxicity. 10,11 
The presence of intermediate aggregation 



www.landesbioscience.com 



Prion 



339 




Figure 1. Molecular model of an A(3 40 fibrillar octamer interacting with a POPE bilayer membrane 
(lateral view). The N-terminal p-strands are shown in red, the core C-terminal (3-strands are blue, 
and the turn regions are green. The hydrophilic lipid headgroups (lime) and the lipid tails (gray) 
are illustrated by spheres. 



products (A(3 oligomers) in the extracel- 
lular or intracellular environment, rather 
than that of fully formed fibrils, seems to 
alter the membrane structure and func- 
tions and is related to the dysregulation of 
cellular ionic homeostasis. 12 

Although the details of the molecu- 
lar mechanism of membrane permeation 
are unknown, amyloid oligomers have 
been shown to interact with lipid mem- 
branes and cause nonspecific ion leak- 
age, 13 " 15 or may transform into annular 
protofibrils that can form ion channel-like 
structures. 16 " 18 

Finding detailed molecular informa- 
tion about the structural organization of 
amyloid fibrils, fibrillar and pre-fibrillar 
oligomers and annular amyloid chan- 
nels remains a major yet challenging goal 
of both experimental and computational 
modeling studies of A(3 aggregates. In a 
recent study, 19 we used molecular model- 
ing in conjunction with atomistic molecu- 
lar dynamics to investigate in detail the 
early steps in the mechanisms of interac- 
tion between preformed A(3 40 fibrillar 
oligomers (generated using protofilament- 
like structures 20,21 ) and a model bilayer 
membrane consisting of 98 POPE lipids. 
The initial protofilament structures are 
built based on experimental and com- 
putational studies of A(3 protofilaments 



formed in solution, 20,22 which have been 
simulated and studied previously in detail 
both under physiological conditions 20 and 
at elevated temperatures. 21 Consistent with 
previous studies, the first eight residues 
from the N-terminus of the A(3 peptides 
are not included in the protofilament mod- 
els because they were shown to be signifi- 
cantly more disordered than residues in 
the 9—40 region. 19 " 21,23 " 25 The main focus 
of our new work 19 is on the atomistically 
detailed analysis of the role of different 
molecular and structural elements in the 
stability and conformational dynamics of 
A(3 protofilaments near lipid bilayers, and 
on the implications on the mechanism of 
possible insertion of fibrillar oligomers and 
protofilaments into the membrane. 

Molecular Interactions 
with POPE Lipids 

In reference 19, we used MD simulations 
to probe with atomistic-level detail the 
mechanism of interaction between molec- 
ular models of protofilaments of A(3 40 pep- 
tides and a POPE lipid bilayer in explicit 
water. The POPE lipids were chosen simi- 
larly to previous studies, 26,27 and because 
the smaller headgroups of POPE lipids 
are expected to promote a more favor- 
able interaction with amyloid peptides. 



Note however, that we have later used 
POPC lipids as well (see below) which 
lead to generally less favorable interactions 
with inserted peptides, being thus more 
stringent tests of the stability of amyloid 
channels. Here, we built protofilament 
systems consisting of eight A|3 peptides 
and explicit water molecules, with a total 
number of 42,000-57,000 atoms, with 
different initial conformations, and simu- 
lated them for 150 ns each in the NPT 
ensemble (i.e., constant number of mol- 
ecules, pressure and temperature), using 
periodic boundary conditions, at a con- 
stant pressure of 1 atm and at the temper- 
ature of 310 K, close to the physiological 
values, similar to previous studies. 19 " 21,28 
The simulation range was long enough 
to observe the first detailed steps of the 
intermolecular interactions. After the A(3 
protofilaments lost their initially ordered 
(3-sheet-rich structures, they reached 
quasi-equilibrated states, preserving only 
a fraction of their characteristic |3-sheet 
content. The change of structure near the 
membrane was significantly greater than 
the change observed for corresponding 
simulations in bulk water, 20,21 due to the 
presence of the specific lipid environment. 

Analysis of both the change in the 
number of peptide-related hydrogen bonds 
(NHBs) and of the electrostatic interac- 
tions between the peptides in the proto- 
filament and the membrane demonstrated 
that the changes observed in structure 
of the protofilaments are due to specific 
interactions with the lipids. 19 We identified 
the electrostatic attraction between the 
charged residues and the lipid headgroups 
acting as the main driving force that leads 
to a further penetration of the C-termini 
into the hydrophobic lipid tails region, as 
illustrated in Figure 1. At the same time, 
we observed perturbations in the mem- 
brane as an effect of the proximity of the 
protofilaments, as the lipid headgroups 
reposition themselves around the charged 
residues, as shown in Figure 2. This fig- 
ure shows that, along the MD trajectory 
there is an area consistently avoided by the 
center of mass (COM) of the phosphorus 
atoms, which coincides with the COM 
positions of the atoms in the C-termini. 
Also, the tails in the upper membrane 
leaflet reorient themselves closer to the 
surface and the C-termini, which leads 
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Figure 2. Repositioning of lipid headgroups due to interactions with the A(3 40 protofilament seg- 
ment facilitates the access of the hydrophobic C-terminal (3-strands into the lipid tail region (top 
view; see also Fig. 1). Here, four top-view periodic images of the simulation box are shown for a 
typical trajectory. The dots represent the positions sampled by the center of mass (COM) of the 
lipid P atoms (pink) and the COM of the C-termini (blue) during the 45-150 ns trajectory segment 
when the C-termini are associated with the lipid head groups. A representative A(3 9 _ 40 octamer 
structure is shown in only one periodic image. 




to a local membrane-thinning effect. 1 '' 
Our study demonstrates that in the prox- 
imity of lipid membranes the protofila- 
ments adopt structures that have a lower 
content of (3-sheet compared with the 
ones formed in bulk water. Although the 
timescales accessible to MD studies can- 
not provide us with a full insight into the 
complete insertion process of aggregated 
A(3 peptide structures into the membrane, 
the early interactions observed show that 
the structure formed in solution cannot 
be transferred to the membrane environ- 
ment without a significant loss of (3-sheet 
content. This concept has direct potential 
implications in the understanding and 
further modeling of A(3 ion channel-like 
structures formed in membranes. 

In our MD studies, after the initial 
protofilament association with the lipid 
bilayer (i.e., at about 45 ns of NPT dynam- 
ics), the content of (3-sheet structure was 
stable for the rest of the simulation and we 
observed little helical structure formation. 
These results suggest that, if the chan- 
nels are formed through direct membrane 
incorporation of (3-sheet-rich fibrillar 
oligomers, these would need to undergo 
further conformational changes before or 
after being incorporated into lipid bilayers, 
and may need to be associated with other 
similar units in order to form (3 -turn- (3 
ion channel-like structures. 

Note also that experimental studies 
showed that it is the prefibrillar, globular- 
like oligomers rather than the fibrillar ones 
that might associate into annular struc- 
tures, 28 but their exact secondary structure 
is still unkown. Another possibility is that 
the predominantly helical membranar 
A(3 peptides may be associated with other 
similarly-structured A(3 monomers before 
leaving the membrane, as there have been 
studies 29 supporting the idea of dimeriza- 
tion in the transmembranar domain of 
APP. On the other hand, experimental 
studies analyzing the structure of A(3 pep- 
tides in membranes, amyloid channels or 
pore-like structures showed that they con- 
tained a high percentage of (3 -sheet, 30 " 32 
suggesting that, if it the aggregation starts 
from a helical structure, the N-terminus 
may adopt (3-sheet conformations, while 
the C-terminus has been shown to have a 
high percentage of helix when in a lipid 
environment. 33 



Polymorphism of A|3 Peptide 
Channels 

The ability of amyloid oligomers to per- 
turb cellular membranes, whether it is 
by thinning the membrane or through 
ion channel formation, is related to their 
cytotoxicity. Our study of A(3 4(J interact- 
ing with a POPE bilayer focused on the 
mechanism of interaction between fibril- 
lar oligomers and lipid bilayer models. 
Understanding these interactions opens 
new questions on the molecular pathways 
of aggregation of A(3 peptides into annu- 
lar channel-like membrane structures. 
Membrane disruption by non-fibrillar 
(globular) amyloid proteins is followed 
by an increase in membrane conductance 
either through specific 18 or non-specific 
ion leakage. 13 The mechanism of globu- 
lar oligomers inducing cell dysfunction 
due to membrane permeabilization has 
been documented in several studies, 14,34 
but their atomistic structure in membrane 
is still unknown. Experimental evidence 
(e.g., X-ray fiber diffraction and solid-state 
NMR) has characterized the A(3 fibrilar 
aggregates in solution as being definitely 
|3-sheet rich, whereas their conforma- 
tions in lipid bilayers is largely unclear. 



Previous studies predicted A(3 ion chan- 
nel conformations based on the peptide 
sequence and on how it relates to other 
patterns observed in different membrane- 
bound and channel-forming peptides 
and proteins. 35 More recently, electron 
microscopy images revealed "pore-like" 
annular structures for amyloidogenic 
protofibrils (i.e., mutants of A(3 4 and 
of ct-synuclein). 17 However, these struc- 
tures were not formed in membranes and 
were not associated within bilayers, their 
atomistic details in lipids being uncertain. 
Using a series of experimental techniques 
(AFM, CD, gel electrophoresis and elec- 
trophysiological recordings), reference 
36 showed that amyloid proteins includ- 
ing A(3 40 , a-synuclein and amylin can 
undergo structural conformation changes 
and associate into ion-channel-like struc- 
tures comprising typically of a few sub- 
units. It has also been suggested that an 
a-helical form can be involved in a process 
of membrane poration. 37,38 

Based on similar experimental stud- 
ies 39 " 42 Nussinov et al. pioneered the atom- 
istic modeling of these structures. 43 " 47 
These models comprise of 10-24 A(3 4 
or A(3 4 , monomers, each with a (3-turn-(3 
structure — the main tertiary structure 
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Figure 3. Possible models of A(3 dodecameric ion channels. The 12 peptides are colored by resi- 
due type: blue, positively charged; red, negative; green, polar; white, nonpolar. Top (A, C, E and G) 
and lateral (B, D, F and H) views of several models (see text). 



motif that characterizes the peptide as 
part of mature fibrils in formed in water. 
MD refinement of these models shows 
that the initially continuous fibrillar 



structures break in about the same number 
of oligomeric subunits as reported by AFM 
imaging experiments, the final channel- 
like structure having similar dimensions 



compared with the AFM images. The A(3 
N-termini maintain essentially their initial 
(3-sheet structures, whereas the C-termini 
become more distorted. These studies, 
together with others looking at structures 
that aggregated A(3 peptides can adopt in 
lipids or different hydrophobic environ- 
ments (e.g., some that consist of two helical 
regions of residues 8-25, and 28-38 33 ' 48 ' 51 ) 
led us to infer that the (3 -turn- (3 might not 
be necessarily the only structural motif 
adopted by A(3 peptides when embedded 
in lipid bilayers. Some of our recent mod- 
els shown in Figure 3 consist of 12 mono- 
mers with a combined (3-strand in the 
N-termini region and helical C-termini. 
Here, 12-mer channels are shown for illus- 
trative purposes, and to keep the size of the 
overall atomistic model small. However, 
similar atomistic models may be easily 
extended to larger (e.g., 16-mer, 20-mer, 
etc.) amyloid-channel-like structures. 
Here, we have also modeled two structures 
containing only ct-helices and no signifi- 
cant (3-strand content. We note that, while 
these all-atom models were built using a 
standard force fields and have been mini- 
mized and equilibrated using explicit water 
molecules (i.e., similarly to the methods 
described in ref. 19) they have not been 
yet fully tested and are only suggestions of 
such possible structures. 

For example, to obtain the first model 
(shown in Fig. 3A and B), we merged 
structures of an A(3 monomer in a two- 
stranded fibril of Tycko's model 22 (which 
has thus originally a (3 -turn- (3 pattern) 
with the helical structure of A(3 40 suggested 
in reference 52. The result, represented in 
Figure 3A, is a transmembrane-spanning 
channel-like structure with the N-termini 
containing charged residues (Dl, E3, R5, 
H6, D7, Ell, H13, H14 and K16) in the 
interior, as (3-strands, and the more hydro- 
phobic C-termini modeled as helices 
toward the exterior of the channel, inter- 
acting with the hydrophobic lipid tails. 
Note that the charged residues are struc- 
turally located such that the first two nega- 
tively-charged residues (Dl and E3) would 
interact with the positively-charged amino 
groups of the lipid headgroups and thus 
stabilize (anchor) the protein segments in 
the lower membrane leaflet. The same kind 
of rationalization goes for the top part of 
the channel, with the charged residues E22 
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and D23 performing the same role. Here, 
we used POPC lipids, which are found 
abundantly in the membrane of neuronal 
cells, have the same tails as POPE, 19 but 
the ethanolamine is replaced by choline in 
their headgroups. Our lipid bilayer has a 
thickness of about 38 A (i.e., between the 
planes containing phosphorus atoms in the 
two membrane leaflets) . 

In Figure 3C, E and G (top views) 
and 3B, D and F (lateral views) are 
shown other possible structural mod- 
els of transmembrane amyloid channels. 
In Figure 3C and D is shown an A(3 4 
channel modeled as a helical structure, 
starting from the monomer conformation 
in aqueous solution of fluorinated alco- 
hols 49 (Protein Data Bank code II YT). 
Each monomer has two helices, one in 
the region of residues 8-25 and at residues 
28-38. Thus the possibility that charged 
residues Dl, E3 and R5 may act as "elec- 
trostatic anchors" to keep the peptides in 
the lipid bilayer. 

In Figure 3E and F is shown an A|3 4 
channel that is essentially helical, modeled 
starting from the PDB structure 1BA4. 52 
The N-terminus is mainly unstructured, 
with the C-terminus forming a long helix 
in the region of residues 20-39. The 
hydrophobic C-terminus can interact with 
the lipid tails, and the charged residues 
can also anchor the structure in the upper 
membrane leaflet. 

Note that these first three cases (i.e., 
shown in Fig. 3A— F) inherit three suc- 
cessive (G-XXX-G) motifs from the 
transmembrane region of the amyloid 
precursor protein. This type of motif is 
known to promote dimerization of poly- 
peptides via C -H- • -O hydrogen bonds 
between two segments in a membrane 
environment 29,53 — a phenomenon that 
might occur when the helical parts of the 
peptides in the amyloid channel interact 
inside the lipid membrane. This factor 
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the molecular structures of membrane- 
absorbed fibrillar oligomers and amyloid 
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